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Abstract Recent progress is summarized on the determination of the density distri- 
butions of stars and dark matter, stellar kinematics, and stellar population properties 
in the extended, low surface brightness halo regions of elliptical galaxies. With in- 
tegral field absorption spectroscopy and with planetary nebulae as tracers, velocity 
dispersion and rotation profiles have been followed to ~ 4 and ^5 — 8 effective 
radii, respectively, and in M87 to the outer edge at ~ 150 kpc. The results are gen- 
erally consistent with the known dichotomy of elliptical galaxy types, but some 
galaxies show more complex rotation profiles in their halos and there is a higher 
incidence of misalignments, indicating triaxiality. Dynamical models have shown a 
range of slopes for the total mass profiles, and that the inner dark matter densities 
in ellipticals are higher than in spiral galaxies, indicating earlier assembly redshifts. 
Analysis of the hot X-ray emitting gas in X-ray bright ellipticals and comparison 
with dynamical mass determinations indicates that non-thermal components to the 
pressure may be important in the inner ~ 10 kpc, and that the properties of these sys- 
tems are closely related to their group environments. First results on the outer halo 
stellar population properties do not yet give a clear picture. In the halo of one bright 
galaxy, lower [a/Fe] abundances indicate longer star formation histories pointing 
towards late accretion of the halo. This is consistent with independent evidence for 
on-going accretion, and suggests a connection to the observed size evolution of el- 
liptical galaxies with redshift. 



1 Introduction 

There are indications that the outer halos of bright elliptical galaxies formed later 
and through different processes than their inner parts. It is clearly important to estab- 
lish whether this is true for some, most, or all ellipticals. Signatures of the formation 
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processes may be preserved longer in the outer halos than in the well-mixed inner 
parts, because of the longer dynamical time-scales at large radii. Owing to the low 
surface brightness in these halos, the outer kinematics and stellar population prop- 
erties can be best studied in local galaxies. Combining the results with predictions 
from simulations and with high-redshift studies of ellipticals one may expect to 
obtain a global picture of halo formation in elliptical galaxies. The following sec- 
tions will discuss some recent results on the stellar distribution and kinematics in 
the outer halos (§2), the dark matter distribution and dynamics (§3), and finally the 
stellar populations and assembly history of nearby elliptical galaxy halos together 
with related results from high-redshift studies and simulations (§4). 

The subject of the outer halos of ellipticals is one that has benefitted enormously 
from the vision of Ken Freeman, whose scientific work we celebrate with this con- 
ference. Ken was one of the first to realize the potential of using planetary nebulae 
(PNe) as kinematic tracers in regions where the surface brightness is too low for 
traditional absorption line spectroscopy. In 1987, in the course of a project with 
X. Hui, H. Ford and M. Dopita, he measured the first PN radial velocities in the 
halo of Centaurus A, the nearest elliptical galaxy, using a multi-fiber instrument at 
the Anglo-Australian Telescope. Together with later AAT and CTIO mea surements, 
they were able to obtain 433 PN velocities in this galaxy dHui et al Together 
with M. Arnaboldi and others, he extended this work to the fainter PNe in the halos 



of more distant ell ipticals such as NGC 1399, M86, and NGC 1316 dArnaboldi et al 



1994 1 19961 1998), in order to constrain the kinematics and dark matter mass distri- 



bution in these galaxies. 

The potential shown by these papers led to the idea of a special-purpose instru- 
ment, the Planetary Nebula Spectrograph (PNS), which is based on the principle 
of using counter dispersed imaging to simu ltaneously find PNe and measure their 
radial velocities dDouglas et al.ll 19971 120021) . This instrument is currently mounted 
at the William Herschel Telescope in La P alma. As two h ighlights, some 2600 PN 



velocities were measured in nearby M31 dMerrett et al.l 120061). and within a still 
ongoing project ~ 100-200 P N velocities have so far been obtained for each of a 
sample of some 15 ellipticals dCoccato et al 1l20Q9h . The 1996 paper on M86 also re- 
ported three discrepant PN velocities, ascribed to the intracluster stellar population 
in Virgo. Since then, Ken has enthousiastically participated in using PNe to stud y 
the dynamics of intracluster stars in nearby galaxy clusters ( Freeman et al. 2000l) : 
see also M. Arnaboldi's paper in this volume. For the firs t multi-slit imaging obser- 
vations of PNe in the Coma cluster ( IGerhard et alj|2005l) he went observing at the 
Subaru telescope. Ken is a coauthor of some 30 refereed papers relat ed to the use of 



PNe f or kinematic measuremen ts, including also the Galactic bulge iBeaulieu et al 



i-Tse t or kinematic measurements, including also 
2000l) and, again, Centaurus A dPeng et al.ll2004 
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2 Density distribution and kinematics of outer halo stars 



Accur ate brightness profiles over large radius ranges were constructed bv lKormendv et al 
d2009h "or Virgo cluster elliptical and spheroidal galaxies. The brightest elliptical 
galaxies in this sample have profiles characterized by large Sersic indices n ~ 4 — 12, 
large effective radii R e > 10 kpc at faint surface brightnesses, [i e y — 23 — 25, and 
central cores (depressions relative to the Sersic profile). Fainter ellipticals have 
Sersic indices n ~2-4, R e 's of a less than a few kpc, ji e y — 19 — 22, and central 
cusps (extra light relative to the Sersic profile). While the dichotomy in the central 
profiles of both classes is clear and is their primary distinction, there are some tran- 
sition cases in the large-radius profiles. The important point for this talk is that the 
(bright) core ellipticals have the most extended luminous halos, an d that within each 
of the two classes n (extent) does not correlate with luminosity dCaon et al 



Kormendv et al 



es n ( e: 
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1993 



Ou ter halo kinematics based on PN velocities were analyzed by ICoccato et al 



br a sample of 15 early-type galaxies. These data typically reach to ~ 5R e 
and in some galaxies to 8R e . The derived profiles of PN number distribution and 
kinematics agree with surface brightness and absorption line kinematics where the 
data sets overlap, within their respective errors. The radial profiles of mean rms 
velocity v,^ for this sample fall within two groups, with part of the galaxies having 
slowly decreasing v n ns(^) an d the remainder having steeply falling v^^R). There 
is a large overlap of these groups with the c ore and cusp ellipticals, respectively 
dKormendv et alj|2009t ISopkins et al.ll2009aUbb . The outer halos of ellipticals tend 
to be more rotationally do minated than the centra l parts. The slow-rotator versus 
fast-rotator classification of Emsellem et al.l (12007b based on inner Sauron data is 
largely preserved in the halos; however, some more complex profiles of specific 
angular momentum parameter Xr are seen in the halos. Slow rotators typically reach 
Xr ~ 0.1 in their outer parts. Lastly, twists and misalignments in the velocity fields, 
commonly interprete d as signs of triaxiality, seem to be more frequent at large radii 
dCoccato et alj|2009l). 



Proctor et al.l (120091) have pioneered a technique to measure absorption line kine- 



matics for the integrated halo stellar population using data from multi-object slitlet 
spectroscopy, and to construct 2D velocity fields to 3R e from these data. They find a 
variety of stellar rotation profiles beyond lR e , constant, decreasing, and increasing, 
confirming that the fast-slow rotator classification may be more complicated in the 
halos than within R e . Halo kinematics obtained with integral fi eld units (IFUs) such 
as Sauron and VIRUS -P rea ch somewhat further, to ~ 4R e ( Murphy & Gebhardl 
2007; IWeijmans et al.l 120091) . In these studies, the signal from many lenslets or 



fibers is co-added to increase the signal-to-noise. Because of the limited field-of- 
view, several IFU fields may need to be exposed. For the galaxies where kinematics 
is available both from these absorption line techniques and from PNe, the veloc- 
ity dispersion profiles agree within the errors. The integrated light techniques yield 
measurements also of the line profile shape parameters, giving valuable constraints 
on the anisotropy of the halo orbit distribution. These are harder to obtain from PNe 
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unless large samples can be obtained. Since PNe can be found out to larger radii, 
the best dynamical constraints will come from combining both techniques. 

Th e extreme outer halo o f the Virgo-centric giant elliptical galaxy M87 was stud- 
ied bv lDohertv et al.l (12009b . They used small samples of PNe in fields at 60 and 150 
kpc to measure the velocity dispersion, after removing Virgo intracluster stars. PNe 
centered around the systemic velocity of M87 were found only out to 150 kpc, but 
not further, and the velocity dispersion profile was found to decline steeply towards 
the outer edge. Given that the X-ray-determined mass profile is rising steeply at 
these radii and is dominat ed by the clus t er, thi s is dynamically possible only if the 
halo of M87 is truncated. McNeil et al. ( 2010b studied the central bright galaxy in 
the Fornax cluster, NGC 1399. Again using PNe, they found that a component of 
low-velocity PNe (by ~ 800 kms ) is superposed on the main galaxy population. 
When this component is removed, the velocity dispersion profile of NGC 1399 is 
flat at ~ 200km s _1 from PNe and red globular clusters out to 80 kpc. Both studies 
show that the outermost halos of elliptical galaxies are in interaction with the envi- 
ronment, and that discrete tracers are important to disentangle the different (mixed 
and unmixed) components in the velocity distribution. 



3 Dark matter and dynamics 



The dark matter distribution in elliptical galaxies can be studied with a variety of 
techniques, including strong and weak lensing, X-ray emitting hot gas, and dynam- 
ics. The lensing results are discussed in K. Kuijken's talk in this volume; in brief, 
while with weak lensing data the dark halos of elliptical galaxies can be followed 
to hundreds of kpc, the strong lensing analyses give accurate mass measurements 
within the Einstein radius, and together with central velocity dispersions constrain 
the slopes of the mass density profiles to an average value of very nearly —2 (isother- 
mal), with small scatter. 

Early dynamical analysis of stellar velocity dispersions and line profile shape 
data with non-parametric spherical models first indicated that the inner circular ve- 
locity curves (CVC) of elliptical galaxies are flat to ~ 10%, out to approximately 



1-2 R e dKronawitter et alfeOOOt IGerhard et ap200lb . Recent work by lThomas et al 



d2007h on a sample of ellipticals in the Coma cluster with slightly more extended 
data (1-3 R e ) has shown more varied CVC slopes, with some falling, some flat, and 
some rising. Using axisymme tric Schwarzschild mo dels to constrain the dark matter 
halos in the Coma ellipticals, iThomas et al. I d2009l) found that the dark matter den- 
sities in these galaxies are on average 7 x higher than in spiral galaxies of the same 
luminosity, and 13 x higher than in spirals of the sa me baryonic m a ss, co nsistent 
with the earlier analysis of round nearby galaxies by iGerhard et al.l (1200 lb . Bary- 
onic contraction is not sufficient to explain the difference, implying that the inner 
halos of ellipticals presumably formed earlier than those of spiral galaxies. If one 
assumes that the halo densiti es measure the density of the Universe at the time of 
collapse (IGerhard et al ] l200lh . the assembly redshift (1 +z a s) for the Coma ellipti- 
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calsis twice that for typical spiral galaxies, and puts their z a 
20091) . with brighter ellipticals assembling later. 



1...3 (IThomas et al 



Several intermed iate-luminosity ellipticals have steeply falling o uter PN velocity 
dispersion profiles ( IRomanowskv et alj|2003[ ICoccato et alj|2009l) . The prototype 
galaxy for this class is NGC 3379, for which some 200 PN velocit ies are available. 
It was thus argued (IRomanowskv et alJl2003b iDouglas et alj|2007l) that this galaxy 
might be 'naked' , i.e., lack a signi f icant dark matter halo. Extensive recent dynam- 



ical modeli ng bvlde Lorenz 



et al 



d2009h . using the made-to-measure particle code 



NMAGIC (Ide Lorenzi et al.l 12007). has shown that a variety of dark matter halos 
are allowed by the data, implying radially an isotropic orbit distribu tions, and that 
the no-dark matter halo is ruled out (see also IWeiimans et al. I l2009h . However, the 
allowed halo density distributions imply that the CVC in this galaxy falls by at least 
25% between the center and 2R e , i.e. they are of significantly lower density than 
would be required for a flat CVC. Similar results were found for two other galaxies 



of thi s class, NGC 4697 and NGC 4494 dde Lorenzi et al.l 120081: iNapolitano et al 
20091) . 



There is a long history of mass determination in X-ray bright elliptical galaxies, 
using the hot g as as a tracer in assumed hydrostatic equilibrium in the gravitational 
potential (e.g. [Nulsen & Bohringer] 1 19951 : Humphrey et al. 1 120061: Fukazawa et al 



120061: Nagino & Matsushitall2009l) . Based on a comparison of the gra vitational po- 
tentials and c ircular velocities inferred from X-rays and from dvnamics JChurazov et al 
J2008L 120 1 Oh argued that the non-therm al sources of pres sure contribute on aver- 
age 20-30% of the thermal gas pressure. iDas et al. (2010) recently devised a non- 
parametric Bayesian technique to determine CVCs and confidence ranges from de- 
projected hot gas pressure and temperature profiles. Comparing with dynamical 
mass determinations, they also find differences indicating that non-thermal pressure 
sources are important in the central 10 kpc of their analyzed six galaxies. Other pos- 
sible interpretations for these differences are violations of hydrostatic equilibrium or 
biases introduced by the assumptions in the dynamical models; this is an important 



Das et al 



Das et al 



issue t hat needs to be clarified. The dark matter mass fractions inferred by 
(201C) are 35 - 80% at 2R e , rising to 80 - 90% at the outermost radii. 
(120 101) also find that for all of their galaxies the outer CVCs are rising, implying 
steeper than isothermal mass profiles. Furthermore, the inferred circular velocities 
and the luminosities of their galaxies are found to correlate with the velocity dis- 
persion of the environment. Both results suggest that the properties of X-ray bright 
ellipticals are closely related to their group environments. 



4 Stellar population and assembly history 

Studies of the stellar populations in the outer halos of early-type galaxies may give 
important constraints on their assembly history. Due to the faint surface brightnesses 
in the outer regions this work is only in its beginning, requiring large telescopes 
and/or special techniques. 
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In a recent study, ICoccato et al 



d2010al) used deep medium resolution spec- 
troscopy with FOCAS at the Subaru telescope to measure Lick line-strength indices 
far into the halo of one of the two Coma b rightest cluster galaxie s (BCGs), NGC 
4889. Combining these with literature data. ICoccato et alj ( 1201 Obi) constructed ra- 
dial profiles of metallicity, [a/Fe] abundance ratio and age, from the center out to 
~ 60 kpc (<~ 4-R e ). These profiles show evidence for different chemical and star 
formation histories for stars inside and outside ~ l.2R e = 18 kpc radius. The in- 
ner regions have a steep metallicity gradient and high [a/Fe] at ~ 2.5 solar value, 
pointing to a rapid formation process. This is consistent with a quasi-monolithic, 
dissipative merger collapse at early redshifts dKobavashill2004l) . followed by one 
or at most a few dry mergers between several such units so tha t the original steep 
metalli city gradient (Chiosi & Carrarol 2002 ; Pipino et al. 2008 ) was only partially 
erased ddiMatteo et alj|2009h . In the halo of NGC 4889, between 18 and 60 kpc, 
the metallicity is near-solar with a shallow gradient, while [a/Fe] shows a strong 
negative gradient, reaching solar values at 60 kpc, and the inferred ages are 9-13 
Gyr. This argues for lat er accretion of st ars from old systems with more extended 
star formation histories ( Matteuccilll994l). presumably through mino r mergers from 
th e galaxy's group env ironment (lAbadi et al hoodlNaab etal]l2009h . 

Rudick et al. (2010) study the colour distribution in the outer halo of M87 and 



the adjacent intracluster light. They find that the M87 halo becomes bluer with ra- 
dius to the outermost radii, and argue that the common colours of the outer halo and 
of some surrounding intracluster tidal features suggest th at the galaxy's envelope 
may have formed fr om similar streams. Foster et al.l (|2009l) used the technique from 



Proctor et al. (2009) to study three galaxies. Their most extended dat aset is for NGC 



2768 w hich they suggest, based on comparison with simulations bv lHopkins et al 



d2009al) . is consistent with a dissipative merger origin. The two other d a tasets have 
too large errors beyond R e to allow firm conclusions. IWeiimans et al. I d2009b find 
from co-added Sauron field data that the metallicity gradients in the two intermedi- 
ate luminosity ellipticals NGC 3379 and NGC 821 remain approximately constant 
to 4R e , with old and subsolar metallicity populations in the halos. Thus outer halo 
stellar population studies to date do not yet c ombine to a clear pictu re. 

The bimodal stellar populations found bv lCoccato et al.l (1201 Obi) in the core and 
outer halo of NGC 4889 are consistent with recent results on the size evolution 



of ea r ly-type galaxies (ET Gs) with redshift (e.g iDaddi et al.l l2005t iTrujillo et al 



2007t ICimatti et al.l 120081) . such that ETGs at z ~ 1 (z ~ 2) have sizes a factor 



of 2 (3 — 5) smaller than ETGs with similar mass toda y dvan der Wei et alj 2008 



van Dokkum et al.ll2008t ISaracco et al.ll2009l) . Recently. Ivan Dokkum et al.l ( 1201 



38: 
10) 



find from stacked rest-frame R-band images that massive ETGs have nearly con- 
stant mass inside 5 kpc with redshift, but increase their envelope mass by a factor 
4 since z = 2, with effective radius evolvin g as R P °c (1 -j-^) - 1 - 3 . Minor merg- 



ers may play an important role in this process (NaabetaL 20091) . and the accreted 



galaxies are likely to have gone through longer star formation histories, leading to 
less a-enriched stellar populations as was found in the halo of NGC 4889. Further 
work on the stellar populations in the outer halos of local elliptical galaxies may 
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become a valuable tool to understand what part of the elliptical galaxy population 
participates in the observed size evolution. 
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